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Abstract. We consider a cubic nonlinear Schrodinger equation with periodic 
potential. In a semiclassical scaling the nonlinear interaction of modulated 
pulses concentrated in one or several Bloch bands is studied. The notion of 
closed mode systems is introduced which allows for the rigorous derivation of 
a finite system of amplitude equations describing the macroscopic dynamics of 
these pulses. 
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1. Introduction and main result 



In this work we study the asymptotic behavior for < e <C 1 of the foUowing 
nonhnear Schrodinger equation (NLS) 

2 

(1.1) ieatu" = -yAw'' + Vr (^) u' + £'^|u'|^w^ x£R'^,t€R, 

governing the dynamics of a wave field u^{t,-) e L'^{M.'^). Here, k G R and the 
potential Vr — Vr{y) € M is assumed to be smooth and periodic with respect to 
some regular lattice T ~ Z'^, generated by a given basis {Ci, . . . , Cd}, S I^'^i J-G- 

(1.2) Vr{y + -f) = Vr{y), Vy e 7 e T = {7 = ^ti 7iC/ G : 7; € Z}. 

It is well known that if k < the solution of (jl.ip in general does not exist for all 
times, i.e. finite-time blow-ups may occur, of. [25] . 

The equation (jl.ip can be seen as a simplified model of the one considered in 
[7]. There the main motivation was to study, from a semiclassical point of view, 
the dynamics of a Bose-Einstein condensate in an optical lattice, described by Vr, 
cf. [HI [9l [19] for more details. To this end a rescaling of the appearing physical 
parameters yields an equation similarly to (jl.ip , but with an additional non-periodic 
confining potential, which we shall neglect in the following. The parameter e <C 1 
then describes the microscopic/macroscopic scale ratio. The main assumption for 
the analysis presented in [7] has been that the initial data u^(0, •) is supposed to be 
of WKB type and in particular it has to be concentrated in a single (isolated) Bloch 
band E({k) E R. These energy bands describe the spectral subspaces corresponding 
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to the periodic Hamiltonian operator 

(1-3) H;^^:=-'^A + Vr{^), 

cf. Section \2A\ below for more details. In the linear case similar WKB approxima- 
tions have been established earlier in [3l I14| , yielding an approximate macroscopic 
description (i.e. on time- and length-scales of order one) of the highly oscillatory 
solution to However, the question concerning a generalization of the results 

in [7], in particular to the case of multiple bands, has been open so far. Here we 
will answer this question for initial data which correspond to a sum of modulated 
plane waves. 

In order to derive an approximate macroscopic description wc shall proceed by 
a two scale expansion method similar to that in [7]. To this end a detailed under- 
standing of the influence of the nonlinearity is crucial. Indeed we will show that 
the solution to p.ip can be approximated (in a suitably scaled Sobolev space) via 

M 

(1.4) u^t,x) ^ ara{t,x)xi„^ (^;fc™) e'C^"— +0(£), 

m— 1 

for M E N, where the set {{km,im) : m = 1, ...,M} is assumed to form a closed 
mode system, see Definition 12.21 of sufficient high order A. As we shall see A will 
depend on the spatial dimension d. The amplitudes Om are the (local-in-timc) 
solutions to the nonlinear system of amplitude equations 

M 

(1.5) idtOm +^'&m ■'^xa.Tn = ^ l^{p,q,r,m) O-pOLq Qr , m=l,...,M. 

p,q.r = l: 
S(Mp,M(if:Mr) = S(Mm) 

The above system describes a so-called four-wave interaction, also known (most 
prominently in laser physics and nonlinear optics) as four-wave mixing, cf. [Tl[51ll5|. 
By (jl.4p we allow for nonlinear interactions within the same band but also consider 
interactions of different bands. In particular, energy or mass transfer between 
different bands is expected due to the presence of the nonlinearity on the right 
hand side of (|1.5p . To our knowledge this phenomenon has never been studied 
rigorously in the context of Schrodinger type equations. 

One should note that the concept of wave mixing is strongly linked to plane waves 
as considered above. Indeed, if one allows for more general phases (which is possible 
in the case of a single pulse [7]), a rigorous understanding, even in much simpler 
cases, is lacking and in particular one can not expect the nonlinear interaction to 
be maintained on macroscopic time-scales in general. On the other hand we could, 
without any problems, allow for (smooth) higher order nonlinearities as considered 
in [7] . This however would result into a much more involved resonance-structure for 
the nonlinear interactions and in order to keep our presentation simple we restrict 
ourselves to the cubic case. 

Before going into more details, let us briefly mention the following mathemati- 
cally rigorous works which, besides [7], are most closely related to ours: In the 
same equation as (|l.ip is considered but in a slightly different scaling. Similarly, a 
nonlinear Schrodinger type model is derived in [5] from a semilinear wave equation 
with periodic coefficients and in |13j from an underlying oscillator chain model. 
Concerning nonlinear wave interactions, there exist several results (mostly three- 
wave mixing) in the context of strictly hyperbolic systems, see, e.g., [T7| . ITS l [20 l [22] . 
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and in the case of microscopically discrete dynamical systems, cf. [10], [TT] and the 
references given therein. 

The paper is organized as follows: In Section [2] wc introduce the basic notions 
needed in the following, i.e. Bloch bands, mode systems and the closure condition. 
Wc also discuss several illustrative examples of wave mixing. Then in Section [3] 
the formal derivation of the approximate solution is given in detail. In Section [4] 
the obtained formal asymptotic description is rigorously established and our main 
theorem is stated for a closed mode system of order A = 2N + 1, with N > d/2. 
In Section O it is shown that this condition can be relaxed though by introducing 
the concept of weak closure. Finally, in Appendix A, we discuss in more detail the 
underlying Hamiltonian structure of the amplitude equations p.Sp . 

2. Mode systems and resonances 

For our work it is essential to study the spectral properties of 

1 

2' 

since they will basically determine the fast degrees of freedom in our model. 



2.1. Bloch's spectral problem. In what follows we denote by Y the centered 
fundamental domain of the lattice F, i.e. 

r d r 1 1 ^ 

(2.1) y:=<^7eR'^: 7 = E^'^'' ^ 

I 1=1 



2' 2 



By 3^ we denote the d-dimensional torus Ry'p, which is obtained also from Y by 
identifying opposite faces. Note that writing H'^{y) then automatically includes 
periodicity conditions. Moreover, Y* denotes the corresponding basic cell of the 
dual lattice F*. Equipped with periodic boundary conditions 3^* is usually called 
the Brillouin zone and hence we shall denote it hy B = y* . 

Next, consider the so-called Bloch eigenvalue problem [?], i.e. the following 
spectral problem on 3^: 

(2.2) Hr{k)xi{y;k)= Ee{k)xe{y;k), keB,eeN, 

where Ei{k) e R and xeiv) = Xeiv] k) denote the £-th eigenvalue and eigenvector 
of the shifted Hamiltonian operator 

Hr{k) := c-"' yHp,,&'' y = 1 (-iV, + kf + Vr (y) . 

Let us recall some well known facts for this eigenvalue problem, cf. [2l [3l [23l [26] . 
Since Vp is smooth and periodic, we get that, for every fixed k € B, Hr{k) is self- 
adjoint on L'^iy) with domain H^{y) and compact resolvent. Hence the spectrum 
of HY{k) is given by 

spec(i?r(fc)) = {Eeik) ; ^ e N} C K. 

One can order the eigenvalues Ee{k) according to their magnitude and multiplicity 
such that 

Ei{k) < ... < E^{k) < Ei+iik) < ... 

Moreover every Ee{k) is periodic w.r.t. F* and it holds that Ei{k) = Ei{—k). The 
set {Ei{k); fc G S} is called the ^th energy hand. The associated eigenfunctions. 
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the Block functions, Xf (y; k) form (for every fixed fc e S) an orthonormal basis in 
L'^{y). We choose the usual normaUzation such that 

(Xfi(fc),X^2(fc))L2(j;) = J Xl^{y\k)xe2{y;k)dy ^ Se^j.^, £i, £2 e N. 

Concerning the dependence on k e B, it has been shown, cf. [U; that for any 
£ G N there exists a closed subset U C B such that Ei{k) is analytic in O := B\U. 
Similarly, the eigenfunctions Xf sxc found to be analytic and periodic in k, for all 
k G O. Moreover it holds that 

E,^,ik) < E,{k) < Ei+iik), ykeo. 

If this condition is satisfied for all /c G S then Eg^k) is said to be an isolated Block 
band. Finally we remark that 

measZ^ = meas{A: e B \ Ei,{k) = Ee^{k), £1 ^ £2} = 0. 

In this set of measure zero one encounters so-called band crossings. The elements 
of this set are characterized by the fact that Ei{k) is only Lipschitz continuous and 
hence the group velocity 1? := VkEe{k) does not exist. 

2.2. Resonances and closed mode systems. Our goal is to derive an approxi- 
mate description of our model for s <^ I. To this end we shall first introduce several 
definitions needed to do so. 

Definition 2.1. For k e B and £ e N we call ^ = {k,£) a mode and M -.^ B x N 
the set of all modes. 

(i) The grapk of all modes is given by 

g = {{k,EfXk)) : {k,£) eM} cBxR. 

(ii) Given a finite set S ~ {/ii, . . . , fiM : il/ e N} C of modes, we call the 
(ordered) mode system of size A G N generated by S. 

(iii) We further introduce T, : ^ B xRhy 

E(m,„,,...,m„J:= ('^*(-l)^+ifc,„,,^(-l)^+i£;,„Jfc™j') , 

\A=1 A=l / 

where ^* denotes summation modulo F*, and we write tj^^'' :~ S(5^) for the 
corresponding graphs. 

In the following the mapping S will describe the possible nonlinear interaction 
of modes (in every order of e). Note that q'^'' C G and moreover, for any A G N, 
it holds g^!"^ c g^J'^^l Since we are dealing with a cubic nonlinearity. we shall see 
that indeed A G N takes only odd values. 

Definition 2.2. Given a finite set of modes S and a subset T of A^. 

(i) An clement {fi„in ■ ■ ■ iMtoa) ^ called resonant of order A to T, if 

(ii) We say that S is closed of order A, if the group velocity VtEiik) exists for all 
fi = {k,£) £ S and 

(2.3) 4^)n(g\4^)' 
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We infer from the above definition that a single mode /i is always resonant of 
order A = 1 to itself. Note that for any A > 2 however, we obtain 

E : {^i,ni , • ■ • , MmA ) ^ cr ■■= {k,E) e B xR, 

which does not necessarily imply a G G^"^ . As for the condition (|2.3p . it is equivalent 
to saying that for all (/Xmi , • ■ • , MniA ) ^ holds: Either S(^„ij , . . . , G 

{B xR)\ Q, or else there exists a e 5 such that . . . , fJ-mA) — ^(a^) ^ ^.s^^- 

The latter obviously means that (/^mi, • ■ • iMitia) is resonant of order A to S. We 
illustrate these concepts in the examples in Section 12.41 

The condition that V kEi{k) has to exist implies that we can not deal with mode 
systems including band crossings, i.e. which include a = (fc, £) such that k <^U. It 
is known that the higher the dimension d and the higher the band index I, the more 
likely one encounters such crossings. Thus, in terms of practical use one can expect 
our analysis to be restricted to cases where only a few bands with low energies are 
taken into account. 

Remark 2.3. Also note that condition (|2.3p is not equivalent to Q^^'^ r\Q C G^s\ 
since due to multiple eigenvalues Ei{k), we may have modes fx € S and jl ^ S with 

2.3. The system of amplitude equations. In Section[3]we shall derive a system 
of amplitude equations describing the macroscopic dynamics of our nonlinear wave 
interactions. Before going into the details of the derivation let us first discuss the 
general structure of the obtained system (see also Appendix A) . In order to keep our 
presentation simple, we shall from now on consider only the case of non- degenerated 
eigenvalues Ei(k). 

Let iS be a given finite set of modes {^i, . . . , ^^^[} which is closed of sufficient high 
order A. Then the equations obtained for the modulating amplitudes am{t,x) G C, 
for m € {1, . . . , M}, are as in (|1.5|) . i.e. 

M 

idtani i'i^m ' ^x^m — ^ ^ ^{p.q.r,m) ^r- 

S(Mp ^^g ) = S(m7ti ) 

Here 'dm € M is the group velocity corresponding to a given mode — {krm ^m) ^ 
iS, i.e. 'dm ■= ^kEi^^{km), and we further denote by 

jy 

the effective coupling constant Ki^p ,j j, m)7 which in general is complex valued. 

As in the nonlinearity in the amplitude equations is again cubic. It takes 

into account the sum over all p,q,r — 1, . . . , M, for which the resonance condition 
of order A = 3 holds. Explicitly, this is equivalent to the following two conditions 

(2.4) kp ~ kq + kr^ km, Eg^ {kp) - Ei^ (fcj + Ei^ (kr) = Ei,^ (km), 

where the first equation has to be understood as a summation modulo T* . The 
equations (|2.4p describe a so-called four-wave interaction^ i.e. three modes ^p, /ig, /i^. 
being in resonance with a fourth one S S. Clearly, the nonlinear structure on 
the right-hand side of (jl.Sp is such that energy transfer between different modes 
can occur. In other words, even if initially some of the amplitudes am(0, •) are zero, 
they will not remain so in general for times |t| 7^ (sec also Subsection l2.4.1l below). 
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Concerning existence and uniqueness of solutions to the above given amplitude 
equations, we only need a local-in-time result. 

Lemma 2.4. For any initial data (ai(0, •),..., aM(0, •)) G H^(M.'^)'^', with S > 
d/2, the system (jl.Sp admits a unique solution 

(ai, . . . , om) e C"([0, T); H^{W)f' n ((0, T); 

up to some (finite) time T > 0. 



Proof. Since the left hand side of (jl.Sp only generates translations by a constant 
velocity -dm G and thus conserves the L-^(M'^) norm of each am{t, x), the assertion 
of the lemma follows by a standard fixed point argument. □ 



Remark 2.5. Discarding for a moment the nonlinear term in p.5|) we want to 
point out that the remaining linear transport part is much simpler than in the case 
of the full WKB type approximation, as discussed in U (for a single mode only). 
In particular we do not run into any problems due to caustics, since the phase 
functions 

ip„i{t, x) = k„i ■ X - tEt^XKi) 
are globally defined. Also note that in the present work the so-called Berry term 
vanishes, since we do not take into account additional non-periodic potentials, cf. 

mm- 

2.4. Examples. In order to illustrate the abstract concepts defined in the former 
subsections wc shall in the following consider several particular examples of non- 
linear wave interactions appearing in our model. 

2.4.1. Four-wave interaction for three pulses. In this example we shall restrict our- 
selves to a set of modes {/.i,„ = {k„i,im) ■ m = 1,2,3} which is closed of order 
A = 3, at least, and which satisfies the following resonance conditions 

(2.5) h-k2 + ki^k3, Ee,iki)-Ee,{k2) + Ee,{ki) = Ee,{k3). 

(Again the first summation has to be understood modulo T* .) This yields the 
following set of amplitude equations 

idtoi + i'di ■ Vxai = Wi {a)ai + 2k(i^2,i,3) oi 02 as, 

15(02 + i^?2 • Va;a2 = VF2(a)a2 + K{i,2a.^3) «i ^3, 

^ idtas + i-ds ■ V^as = ^3(0)03 -I- 2,1,3) ^2, 

where we shortly denote 



(2.6) 



i',in) l^-j" I 



^^m(^) ■ — ^(m,m,m,m) I'^m I 4" 2 ^ ^ ^{mj,j,' 

We consequently expect the solution of p.ip to be asymptotically described by 

3 

m— 1 



It is easily seen that the nonlinearities W,„(a)a,„ can not transfer energy from one 
band to the other as they are homogeneous in the respective am(t,x) and that all 
of the appearing K[m,j,j,m)j '^■,3 = 1)2,3 are indeed real valued. What is more 
important though is the fact that the above given amplitude system (|2.6p includes 
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an invariant sub-system. Namely, if initially ai(0, •) 0, it remains so for all i G M 
and thus the above given system simplifies to 

idta2 + ii92 ■ ^ = (K(2,2,2,2)|a2p + 2^(2,3,3,2) ksP) a2, 
idta^i + ii?3 • V^az = (^(3,3,3,3)103^ + 2^(2,3,3,2) |a2p) a-^. 

However, such a decoupling does not exist if initially 02(0, •) = 0, since this am- 
plitude will be generated during the course of time by the remaining two others. 
Analogously 03 is generated by interaction of oi and 02- 

More formally we can also infer this fact from the closure condition, since any 
possible combination of ^2 and ^3 via the mapping E (with A = 3) yields either 
^2-1 kz, or any other value (like 2k2 — k^) which is not in Q^^^ by assumption (recall 
that the system of modes ^i,/i2,M3 is assumed to be closed of order A = 3). 
However, if one aims to follow the same argument for, e.g., ki and ^2, the first 
equation in (|2.5[) shows that this sub-system is not closed. 

2.4.2. The case of several pulses within one Block hand. This is a particular situa- 
tion where we keep the band index ^ G N fixed and thus only consider the interaction 
of several pulses within a single Bloch band. Hence the solution to takes the 
asymptotic form 

M 

«^(t,x)~ ^a,„(t,x)x^(^;fc„)ei('=™— ^^('="-)*)/^ + 0(£), ^ S N. 

m— 1 

In particular we can expect this description to be correct in cases where the Bloch 
band Egik) admits only a moderate variation var^ge Eg{k) and is well 

separated from the rest of the spectrum of Hr{k) by a sufficiently large gap, i.e. 

min{|£;^(fc) - E„{k)\ : n G N, 71 7^ ^} > Af > 0. 

It is easy to show then that a four-wave interaction can always be realized within 
such a band. Choose fcmax and k^i^ such that -E^(fcmin) < E^{k) < Ei{kma.x) for all 
k E B. We are looking for a triple (fci, fc2, ^3) satisfying (|2.5p with £j = £, j = 1, 2, 3. 
To this end we note that ~ 2ki — k2, modulo F*, and define 

e(fci, fe) = 2Ee{ki) - Ei{k2) - Ee{2ki~k2). 

Then we have e(fcniax, fcmin) > > e(fcniin, ^max) and by a simple application of the 
intermediate-value theorem, we easily find fci, fc2 with ki ^ k2 and e(fci, fc2) = 0. 

2.4.3. The case of a single pulse decomposed into several bands. This is a second 
particular case, where we expect an asymptotic description for solutions to (jl.l[) 
given by 

L 

u%t, x)^Y. (-' ^0) e'C^-^-^^C^")*)/- + 0{e), fco e M, 

e=i ^ 

where fco is some given wave vector. One should have the following intuition: Given 
an initial plane wave of the form 

<„(x)=/(x)e"=°-/^ 

one decomposes the slowly varying macroscopic amplitude f{x) into a sum of terms, 
each of which is concentrated on a single Bloch band. This is possible due to Bloch's 
theorem, which ensures that L^(R'') = where Tie are the so-called band 
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spaces. Strictly speaking though, one would require countably many terms in the 
decomposition, which we can take into account here. However for any practical 
purposes (and if f{x) is sufficiently smooth and rapidly decaying) only the first 
few Bloch bands need to be taken into account as all higher bands give negligible 
contributions, cf. [16]. We shall not go into more details here on the precise 
definition of Hi etc. but refer to [21 [3l [23l [26] for more details on these classical 
results (see also |16j for a numerical approach). 

Concerning the possible generation of resonant modes, it is clear that in this 
case the first condition in (|2.4p is trivially fulfilled, since we are only dealing with a 
single wave vector fco. Thus, it is merely a question on the precise structure of the 
bands Ei{ko) whether one can indeed expect resonances. 

3. Formal derivation of the approximate solution 

In the following we consider a finite set S of modes which is closed of order 
A = 2N + 1, for some N £ to be determined later. 

For solutions of p.ip we seek an asymptotic two-scale expansion of the following 
form 

(3.1) u%{t,x) := ^e"-Vn(t,x,-,-\ , 

n=0 ^ £ £/ 

where 

(3.2) Vn{,t,x,T,y) -.^ ^ An,a{.t,x,y)Ea {T,y) , 

with E^(T,y) e^'^iy^-r)^ for ^ e 0^^"^+^'^ = S(52"+i). Note that in the most 
simple case, where n = 0, we get that a = {k,Ei{k)), with {k,£) e S, which yields 
Eo-(t, y) = e'^'^'^"^*^'^''^', a simple plane wave. Moreover, if aj G g^^^^~^^\ foj- 
j ~ 1, 2, 3, this yields 

(3-3) EcriEcrjEo-g = Eo-i-o-2+CT3, 

with ai — (72 + 173 G tJ^2(rii+n2+n3)+3) q£ ^^Q^j-gg .^^g j^JgQ ijjrposc that 

3.1. The general strategy. As already said before, we shall only consider the 
case of simple eigenvalues Ei{k). for simplicity. Plugging the ansatz (|3.1|1 . p.2p 
into (jl.ip and expanding in powers of e, we formally obtain 

N 

iedtu% - -ffpor^AT - eK\u%\'^u% = ^ e"X„ + res(u^), 
with the residual 

3Af+l 

(3.4) res(u^)= ^ e"X„. 

n=N+l 

Introducing, for any (general) a = {k, E) E B x IS., the operators 

(3.5) := E - Hr{k), := 19* + ifc • V, + div, Vy, L2 := ^ A. 
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this yields 

(3.6) Xo := (^o^o,.)E, 



and also 



(3.7) 



0-3 



3 = 1,2,3 



where we have used the relation (|3.3p . In general wc get for n = 2, . . . , 3A^ + 1, 
(3.8) 

Xn := ^ (,^o^ri,(T)Ecr + ^ (L^ ^„_i,a-)ECT + ^ L2An-2,a^o 



crey^ ""SSs ""SSs 

— K 



^ ^ ^ ^ ^Tli ,cri^n2,0'2^n3.0'3ECTi— (T2+cr3 ■ 



Here one should note that An,a{t, x, y) = 0, for all n > TV + 1, by assumption. 

Now we shall subsequently construct An^a such that Xn = for n = 0, . . . , A^. 
To this end we have to compare equal coefficients of Eo-. We consequently obtain, 
from p.6p - (|3.8p . equations of the form 

(3.9) LlAn,a = Fn,a, 

where the r.h.s. F„ g. can be determined from the coefficients {A„i g.) ^ri^m+i) for 
771 = 0, l,...,7i — 1. More precisely 

Fn,a '■= ^ L^An-l^a ~ ^2^71-2,(7 

(3.10) + ^ni, 0-1^712,(72 ^"3. '^'3 ■ 

Note that here the summation index has changed in comparison to p.8p . To proceed 
further we need to distinguish two possible cases: 

Case I. On the one hand, for a G Q^^^'^'^^ XQ's^ the closure condition up to order 
2n + 1 implies invertibility of L[^, i.e. (Xq)"^ G Uii{L'^ {y) , H'^ {y)) and wc obtain 

(3.11) AnAt.x.y) = {L^)-^Fn^„{t,x,y). 

The corresponding modes are called non-resonant. 

Case II. On the other hand, if cr G Gg \ then Lg has a nontrivial kernel. In 
order to distinguish this case more prominently from the one above, we shall from 
now on use the notation = cr e which also characterizes the basic resonant 
modes n € S via <^ = S(/i). 

Using the orthogonal projections onto this kernel the necessary and sufficient 
solvability condition for (|3.9p in this case is then given by 

(3.12) F,Fn^,{t,x,y)=0, 
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which yields (recall that dim (ran P,;) = 1, by assumption) 

(3.13) (X., J^n,,)i2(y) = 0, forO^x. eranP,. 
Under this condition, we consequently obtain 

(3.14) An,.;{t,x,y) = an,^{t,x)xdy) + ^n,A*'^^y) 
with 

(3.15) Ai^:^{Ll)-\l-P,)F^^,. 

Note that here a„_,; is still undetermined. However the condition (|3.13[) provides a 
partial differential equation for a„_i^5, the so far undetermined part of the previous 
step. 

Remark 3.1. In the case of non-simple eigenvalues E((k), i.e. dim(ranP^) = R > 
1, we can simply use a smooth orthonormal basis {x<;,r}^=i of ranF^ and generalize 
the above given formulas (|3.12[) . (|3.14p accordingly. 

3.2. Explicit calculations. In the following we shall determine the approximate 
solution in more detail, by following the above described strategy. 

n = 0: We need to solve Xq = and immediately note that in this case. Case I 
above is obsolete. Then, in Case II, equation (|3.6p implies 



L^^o,. = 0, ^eoi'^ 

and thus p.l4p simplifies to 

(3-16) Ao,5 = ao,,;{t,x)xdy), 

with ao,<; still to be determined. 

n = 1: In the next step we have to solve Xi = 0. In Case I, i.e. for a G G^\Gs^\ 
the equations (|3.7p and p. lip imply 

(3.17) ".1-2+.3 

= K ao,<;iao,52ao^<,-3 (Lq) ^ (x?iX?2X?3) > 

fl.i2>53£S^ : 
<r=si -52+^3 

where for the second equality we simply insert p.l6p . We proceed with Case II: To 
this end the solvability condition p.l2p for <; G CJ^^-* allows us to determine the so 
far still unknown ao,^, obtained before. By p.lSp . this yields 



(3.18) / xMFiAt,^,y)'iy^o, 

Jy 

where 

(3.19) Fi^, = L\Ao^, - K Y ^o,,i^,^2^o,< 



^=^1 -52+^3 

From the definition of L\ in p.5p and using the following basic identity 
ixt, (-iVy + k)xe)L2^y^ = VfeS,(fc), 
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a straightforward calculation shows, cf. the appendix of [7], that ()3.18|) can be 
written as 

(3.20) dtao^^ + ■ Vi^ao,.; + ^ i«(;i,?2,<;3,?)'3^o,?iao,?2«o,<;3 = 0, 

'iii^2r';3£y5 = 

''" = ^l-'^2+^3 

where we denote 



''(';i,?2,?3,?) 



(2/)x^2 (2/)x«(y)x<;(y) dy- 



Since any (t g fj^^^ corresponds to a unique m = 1, . . . , M, via = S(/i„i) = 
{km, Ei^{km)), with /im = {km,tm) € 5, we can shortly write 

ao,<;(i,a::) = a„(t, x), X? (y) = X£™ (y; ^m), i?? = "i^m- 

Hence the amplitude equations (|3.20p can be equivalently written in the form (jl.Sp 
used before. 

In summary we have now fully determined the expressions (j3.16p and p.l7p , and 
from p.l4p we finally get that 

where again the coefficients ai^^ are still arbitrary and have to be determined by 
the solvability condition for n = 2. Since this yields an initial value problem for 
ai_^ (sec below) we are free to choose its value at time t = 0. For simplicity we 
shall put ax,^{Q, ■) = 0. 

n > 2: From here we proceed inductively, as described above, by solving X„ = 0. 
The only difference that occurs is that in Case II the coefficients a„_<; do not solve 
a nonlinear initial value problem, but rather a system of linear, inhomogeneous 
transport equations. Indeed, lengthy calculations show that 

dtans + "dci ■ xans + 2iK(f ^^j^ .^^ ;;)a„,fao,?iao,?2 

(3.21) '""^"^^ 



i«:(<;i,?,ft,?)^o,?ia„,5ao,52 + iQn — 0, 



^1 — '^ + ',-2 = '^ 

with source term 

e„ := - P, {L\Ai^ + L^An-i,,) 

f,5i,f2eSs'^ ii,f,'.-2eSs 
^-^1+^2=^ ',-l-^+';2=^ 

+ kP/ 



5i.';2esi;"'?^: 5i,52esi;"'?^: 

cr — i;^ +^2='^ ^1 — CT + ';2=''" 

"l+"2 + "3="^ (2„j+l) 

Again we shall put a„.<;(0, •) = for simplicity, since we are free to choose the initial 
values for p.2ip . Finally, we note that, in order to satisfy = (i.e. in the last 
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step), we do not need to determine the corresponding ajv,? and thus we can impose 
aN,iit, •) = 0. 

4. Justification of the amplitude equations 

In order to obtain our main result. Theorem 14.51 we have to justify the above 
given formal calculations rigorously. In particular we need a nonlinear stability 
result on our approximation. 

4.1. Estimates on the approximate solution and on the residual. In the 

above section we derived an approximate solution which admits an asymptotic 
expansion of the following form 

t X 



u%{t,x)= E ao,?(t,a:;)x? (7) E, 

N 

(4.1) +E^" E {anAt,^)x,{^)+Ai,{t,x,^ 

1) 



n=l ^^^il) 



t X 



N 

"=1 aeG*?"+^>\G(^' 



where the first two lines on the r.h.s. include only resonant modes while the last 
one takes into account the generated non-resonant terms. To this end we required 
the appearing finite set of modes to be closed of order A = 2N+1. 

To proceed further let us, for s G [0, 00), introduce the scaled Sobolcv spaces 

{r e L2(M^) ; < 00} , 

with 

11/11?^= := / il + \ep\'r\m\'dp. 

Here / denotes the usual Fourier transform of / in L^(M''). Note that in the 
following Gagliardo-Nirenberg inequality holds 

(4.2) ys>^3C^>0: || / < / H^, < Coo£-''/'|| / H^. , 

where the factor e^**/^ is easily obtained by scaling. In the following lemma we 
collect the a-priori estimates on needed in order to prove our main result. 

Lemma 4.1. For d,NeN and s e [0, 00) let K = max{0, s + | - 2} and S > 
N + S + ^, or, if d^ 1, S > N + s + 1. Moreover, assume that d'^Vr G i°°(3^), 
for all \a\ < K, and let (ai,...,aA/) S C°([0, T), ^■^'(R''))^^ be a solution of the 
amplitude equations (|1.5p . 

Then, the approximate solution m^, given in (|4.ip . satisfies the following esti- 
mates. For each G (0,T) there exist positive constants Ca, Cb, Cr > 0, such that, 
for £ G (0, 1), |a| < s, and t G [0, T,], it holds 

II {ed.ru'^it, ■) < Ca, II uUt, •) \\h^ < Cb, II res(w^^)(t, •) |1^. < C.E^+i. 
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Proof. Since is given by ()4.ip . in order to prove the estimates of the lemma 
we need to estabhsh the regularity of An. a for n = 0, . . . , A^, cr e G^^^^\ which 
have been calculated formally in Section [3l Under the assumptions (ai, . . . , om) S 
C°([0,T),i/'5(M''))^^ and a"Vr € for \a\ < K, we shall first prove that it 

holds 

(4.3) An.a e C''([0, T), i7^-"(M^ H^+'^{y))). 

Indeed, by (|3.16p . Ao.;;(t, x, y) = aQ t;(t,x)xt;{y), where aq^^ solves p.20p and x<; is 
given by (|2.2p . and wc get immediately (|4.3p for n = (cf. Lemma and Lemma 
lU below). 

From here it follows, on the one hand, that for a e Q^^^ \ g^^^\ by ([SI]) and 
(ISTTll . Fi^„ = io^i,^ e C°{[0,T),H^{R'^,H^{y))) with S = /\ + 2 and hence, by 
(IXTUll . Ai,^ e C"([0,r),ff^-i(R'^,7?^(3^))). On the other hand, for <; € G^^^ we 
obtain bv ' (jTT9)l Fi,<; € C°([0, T), iJ'5-i(K'^, iJ^(3^))), and thus, by (pl^ . Aj^^- € 

C°([0,T),i/^-i(R'',i/^^(3;))). Moreover, ai,,x? e C"([0, T), i?^-i(lK'', ^'(/))), 
since ai^^ solves the linear inhomogcneous transport equation (|3.2ip (for n = 1) 
with initial condition ai.;;(0, •) = 0, coefficients in C'~'{[0,T), (R'^)), and source 
term Gi G C°{[0,T),H'^-^{R'^)). Hence, by (jTU)) . we obtain ^ for n = 1. 
Proceeding inductively, we obtain (|4.3p for all n = 0, . . . , iV via the above given 
steps. 

Having established (|4.3p . we aim to show the first estimate of the lemma. To 
this end we have to guarantee that for each 7i = 0, . . . , TV, ct e G'^^^^\ there exists 
a C > 0, such that 

II {ed^TAn,, (t, ., -) II <C 

holds true for all e E (0, 1), \a\ < s, and t G [0,T*]. Using the Gagliardo-Nirenberg 
inequality (|4.2p . we therefore require 

(4.4) {ed,rAn,a (t, -) e H"\R'') with m > ^, 
for all |a| < s, i.e. 

An,a {t, •, -) e i7f * (R"*) with 771* > S + ^. 

Thus, by ((43)) . we need m* = S - n + s > s + ^, that is, either 5 - ri > + f and 
s > s-|/3|, or S*-/! > and s > s-|/3| + |, for all |/3| < s and aU ti = 0, . . . , A^. It 
turns out that in order to satisfy either of these conditions, the optimal conditions 
we have to impose concerning the regularity of An^a, given in (|4.3p . are S*— ?i > s + | 
and s = K + 2>s+^, for 77. = 0, . . . , TV, which directly yields the assumptions on 
S and K stated in the lemma. 

The second estimate of the lemma then follows immediately by (j4.4p . In order 
to obtain the third estimate, having in mind the definitions (|3.4p and (|3.8p with 
An^a = for 71 > iV + 1 and a^^^ = (cf. the note after (|3.2ip '). it is necessary and 
sufficient to assure additionally that Aajv-i,? e H''{R'^), i.e., ao,^ G i?^+*+i(M'*). 
This yields the additional condition S > N + s + 1, and concludes the proof. □ 

Remark 4.2. Note that in order to derive the above given a-priori estimates for 
UN, the required regularity imposed on is indeed independent of N. 

We shall also need the following result on the linear time-evolution. 
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Lemma 4.3. For e G (0, 1) denote the unitary propagator corresponding to the 
linear Hamiltonian H^^^, defined in ()1.3|) . by 

W{t) := c-'^p-*/^ 

For s G [0, oo) assume 9"Vr G L°°[y) for all\a\ < max{0, s— 2}. Then there exists 
a Ci > 0, such that 

(4.5) II UHt)r Whi < Ci WfWHs for all t e M. 

Proof. Recalling the definition of H^^^., given in p.3|) . we clearly have the basic 
estimate, 

\\f%t)\\^, EE \\U%t)f^\^, = liriL., for all t e M, 
since H^^^ is self-adjoint. Moreover, since [/"^(i) obviously commutes with (its 
generator) H^^^., so does any power of the latter and we therefore obtain 

\\{H^J'/^f%t)\\^, = ||(i/^„r/'.nL., for all s,teR. 

Without loss of generality we assume here that H^^^. > 1, otherwise we may add a 
constant independent of £ £ (0, 1). 

For s = 1 this is nothing but energy conservation. Using the conservation of the 
norm and the condition Vr € L°°{y) we immediately obtain the desired result 
gll) for s = 1 with Ci = (l+4||yr||oo)'/'. 

For general s € N\{1}, we use integration by parts and the assumed regularity 
of Vr to find e-independent constants Ci , C2 > such that 

ciii(^^er)^/'r(t)iL. < iirwii^i < c2\\iH;,,r/'fm)\\^,, 

i.e. we have equivalence of the norms uniformly in e. For general s > the same 
holds true by interpolation. In summary we obtain boundedness of the unitary 
group U^{t) in all and the assertion of the lemma is proved. □ 

4.2. Stability of the approximation. First let us recall the following Moser-type 
lemma, cf. j22l Lemma 8.1], which we shall use in the proof of Theorem 14.51 below. 

Lemma 4.4. Let R > 0, s G [0, 00), and N{z) ~ k\z\'^z with k e K. Then there 
exists a Cs = C's{R, s, d, k) > such that if 

\\iedrfh^<R, V|a|<s anrf||.g|Loo <i?, 

then 

\\m + g)-m)\\H^^<Cs\\g\\Hr 

Proof. The proof can be found in [22] for s G N and follows by interpolation for 
general s £ [0,oo). □ 

With the above results at hand, we are now able to establish our main result, 
which justifies rigorously the validity of the amplitude equations ()1.5|) . describing 
the macroscopic dynamics of M modulated pulses for a closed mode system of order 
A = 27V + 1 (with N depending on d). 

Theorem 4.5. For d e N choose s, 5 e [0, 00) and N eN such that N,s > ^ and 
S > N + s + ^, or, if d 1, S > N + s + 1. Let K ^ max{0, s + | - 2} and 
assume d^Vr G L°°{y) for all \a\ < K. Moreover, let the finite system of modes 
5 = {^1, . . . , //A/ : A/ G N} be closed of order A 2iV + 1. 
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Then for any solution [ai, . . . jOm) G C^{[0,T), {M.'^))^^ of the amplitude 
equations (|1.5p . and any t^ G (0,r), /3 G {^,N], c > 0, there exist an Eq € (0,1) 
and aC > 0, such that for alls € (0,eo); approximate solution u^, constructed 
above, and any exact solution of with 

\\u'{Q,-)-u%_,{Q,-)\\Hs<ce^, 

satisfy 

||u"(t,-)-'a^-i(^r)IL. <C'e^ for all te [0,^*]. 

Remark 4.6. Note, that the approximate solution contains terms up to order 
©(e^), whereas the above given error estimates include only M|^r_l■ The reason for 
this is that for our proof we need to work with but still the obtained error is 
only of the order 0{e^) with (3 < N . We can therefore eventually neglect the last 
term in u^, or, loosely speaking, we can move it to the right hand side of the above 
given estimate. 

Proof. We write the exact solution of (|l.ip in the form 

u'^{t^ x) = u)^{t, x) + e^w^{t, x) 

and denote e'^(i) := || w^(t) Wjjs- Then, clearly, g'^{Q) < c, and we will show that 
there exist C > 0, Eq G (0, 1), such that, for all e £ (0,£o], it holds g'^{t) < C for 
t € [0, i*]- Inserting this ansatz into (jl.ip . written as 

idtu" = -|Au" + -Vr (-) u'+Af(u^), where Af(z) = k \z\'^z, 
2 e \e / 

and applying Duhamcl's formula we get 

w'it)^U'it)w'iO) + e~^ f W{t~T){N{u'j,{T) + e"w'{T))~N{u%{T))) dr 

Jo 

-£-('3+1) / [/"(i-r)res(u^(T))dr, 
Jo 

where res(u^) is defined in equation (j3.4p . Hence, using Lemma 14.11 and Lemma 
14.31 to estimate the residual and the linear semi-group, respectively, we obtain 

Q'{t)< Cic + CiCre''-^'t+ f Ci\\N{u%{t) + e^w' (t)) - N{u'j^{t))\\^^ dr, 

JO 

since g^{0) < c, by assumption. Using N > j3 and e e (0,1), we consequently 
obtain 

g'{t)< Ci{c + Crt,)+e-^ I Ci\\N{u'i,{T)+ePw'-{T))-M{u%{T))\\„^dT, 

Jo 

for t < t^. 

Now, wc set C := C;(c + Crt*)e'-^''^=** and choose a, D > max{c, C}. Then, since 
D > c> g^{0) and g^{t) is continuous, there exists, for every e £ (0, 1), a positive 
time t^jj > 0, such that g^{t) < D for t < if^. 

The Gagliardo-Nirenberg inequality (|4.2p yields, for s > d/2, that 

II e^w'it) 11^^ < e^-''^^C^D for t < fj,. 
Hence, using (3 — d/2 > there exists an Sq € (0, 1), such that 

II e'^w''{t) 11^^ < Ca, for e £ (0,£o] and t < t%. 
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Moreover, by Lemma [47T] we have \\{ed)°'u%{t)\\^^ < Ca for |a| < s, £ G (0,1), 
and t < T. Thus, we can apply Lemma 14.41 (with R = Co) m order to estimate the 
nonhnear term and obtain 

Q'{t) < Ci{c + CrU) + CiCs [ g'{T) dr, for e e (0, Eq] and t<t'h. 

Jo 

Gronwall's lemma then yields 

g'it) < Ci{c+Crt^)e'^'^'-* < C, for e e (0,eo] and t < t^. 

Since C < D, we conclude that the assumptions needed in order to apply Lemma 
14.41 are fulfilled for aU e G (0,eo] and all t < t.^,, that is tf, > f*. Thus the above 
given estimate proves that 

\\u'{t,-)-u%{t,-)\\H^^=Oie^), forte [0,t*]. 

However, since N > (3 and 

\\u%{t,-) - u%_^{t,-)\\^^^ = 0{e''), forte [0,t,], 

we can finally replace u^(t, •) by u^_-^(t, •) in our stability result, which conse- 
quently proves the assertion of the theorem. □ 

We want to stress that our stability result is an advancement when compared to 
the result of [7] , in the sense that our asymptotic estimates no longer suffer from a 
loss of accuracy in powers of e (as has been the case in the cited work). We infer 
in particular, that our nonlinear setting allows for the same kind of stability result 
as one would expect in the linear case. 

Theorem l4 . 5 1 also shows that no blow-up can occur in the exact solution within 
the interval t e [0,r), determined by the existence time of p.5|l . Hence, if the 
system (jl.Sp indeed admits global-in-time solutions we deduce that the solutions 
u'^, starting close to such modulated pulses, exist for arbitrary long times (but 
blow-up may occur after the modulational structure is lost on longer time scales.) 

Finally note that for d = 1 the assumption N > d/2 and the required closure of 
order A = 2N + 1 imply that, at least, = 1 and thus A = 3, whereas for d = 3 
spatial dimensions we need at least N = 2 and hence A = 5. This increase of the 
required A for higher spatial dimensions d can be relaxed though, as we shall show 
in the following section. 

5. The case of higher-order resonances 

In Section [21 we derived the approximate solution under the assumption that 
S satisfies a closure condition of sufficient high order A, i.e. we required A = 2iV-|- 1 
and N > d/2. We shall show now that this closure condition can be significantly 
relaxed if we generalize the ansatz p.ip . (j3.2|) slightly by allowing for a larger set 
of modes Q2N+i- To this end we consider the new ansatz 

N 

(5.1) u'}^{t,x) = Y, e"'An.,a{t,x,y)E^{T,y), 

n=0 creQ-m + i 

where the set of modes ^2A'+i is then defined inductively as follows: Again we start 
from S = {(fci,£i), {hMj^Ai)}, which is now assumed to be only closed of order 
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A = 3 and we set Qi = ■ However, for n > 2 we will allow the sets G2n+i to be 
larger than g^^"+'^\ The reason for this enlargement is motivated by the fact, that 

may be nonempty for some n > 2. Then, the corresponding equation (|3.9p might 
not be solvable. To circumvent this problem, we set 

^ ^(2n-l) , I Pr(2n+1) 

At stage 1 the corresponding An-i.a with a G Q^^"~^^^ take the usual product 
form an-i,a-it, x)xa-iy) with Xa- G kerLg C H^{y) and free coefficients an^i.a- 
Then, in step n the corresponding equations for a e Q^^^^^^^ can be solved as in 
Case II of Subsection 13.11 i.e. we obtain a linear transport equation for the free 
coefficients from the solvability condition for the previous step. 

To make this more precise, we associate with the sequence Gi^Gs, ■ ■ ■ ,G2N+i 
another sequence Gs^Gs, ■ ■ ■ , G2N+1 as follows 

(5.2) ^2n+l = 1^ (^^2rii+l ~ ^2n2 + l + ^2n3 + l^ • 

"1,12. "3=0 

Hence, G2n+i contains all modes of order ©(e") that can be generated by the cubic 
nonlinearity \u\'^u from the already given terms of lower order. 

Definition 5.1. A mode set S = {{ki,£i), (A:m,^m)} is called weakly closed of 
order 2N+1, if there exists a sequence Gi C G3 C ■ ■ ■ C G2N+1 C B xR such that 
the following conditions hold: 

(i) ^2n+i is finite for n = 0, 1, TV; 

(ii) (?i = 4i)=S(5); 

(iii) ^2n+i C fei+i for n = 1, ...,N; 

(iv) if cr G G2n+i n G with n G N, then either cr G G2n-i or a ^ G2n+i\ 

(v) for each cr G G2N+1 n G the group velocity da = \i'kEi{k) exists, 

where ^3 C • ■ ■ C ^2Ar+i is the sequence associated with Gi, G2N+1 according to 
(1521). 

Condition (iv) means that any occurring resonant mode must either occur al- 
ready in an earlier step and hence has a corresponding free coefficient or it appears 
for the first time but it is not yet generated by the nonlinear interaction. This will 
become clearer in the following examples. 

Example 5.2. We first show that a system which is closed of order 2iV+l is 
also weakly closed of order 2N+1. For this we simply set G2n+i = fj^^"'*'^-'. By 
construction, the associated sequence is G2n+i = G2n+i and the conditions (i) to 
(iii) hold immediately. Moreover, (iv) follows since ^2ri+i G = Gi = G^J''^ ■ 

Example 5.3. We now show that a weak closure of order A is in general a 
weaker condition than closure of order A is. For this consider the case S = 
{(/ci,4),(fc2,£2)} such that Gi = G'i^ = {cti,(T2} C G with aj = [k,,Et^{k.j)). 
An easy induction argument gives 

4'"+^) = { (7^+l-J>i + {3-n)a2 ; J = 0, 1, . . . , 2n + 1 } 
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Now assume that G CiGg C {ci, (T2, 3cti— 20-2} for some N* > 2, i.e. we have 

a closure of order 3 but not of order 5. 

We claim that a weak closure of order 2iV+l still holds for any TV such that 
3iV < 2N*. For n = 0, . . . , iV we let 

^2,1+1 = { jcri + (l-j)cr2 ; j = -f3n, ...,«„}, 

with Q!„ ~ [3(n+l)/2] and /3„ — [3n/2], where [•] denotes the integer part. By 
construction of the sequences we thus have Gs ~ Gs 'J {3cri— 2cr2} and 

G2n+i = G2n+i for Ti ~ 2, . . . , N . Hcncc conditions (i) to (iii) of Definition 15.11 are 
fulfilled. Condition (iv) also holds since G2n+i G C Gi ^ {3(Ti— 2(T2}. 

Repeating the construction of Section 13.21 analogously we obtain the following 
result. 

Lemma 5.4. Let S ~ {(fci,£i), . . . , (fcA/,£j\/)} be a set of modes that is closed of 
order A = 3 and weakly closed of order A = 2N+1 for some iV G N. Then an 
approximation 

N 

(5.3) u%{t,x) = J2£" J2 An,At,x,y)'E^ 

n=0 cr&G2n+l 

can be constructed as in Sections \3.1\ and \3.'A where Gi 'Z ■ ■ ■G2N+1 is the sequence 
guaranteed by Definition \5.1\ 

Hence, via the same steps as in Section |4] one can easily obtain a justification of 
the amplitude equations (|1.5p under these relaxed resonance conditions. 

Corollary 5.5. Let d, N Cz N and S he a mode system that is closed of order A = 3 
and weakly closed of order A = 2N + 1. Then, under the same assumptions as before 
the statement of Theorem \4-.5\ holds analogously with the approximate solution given 
by dOl). 

Appendix A. Hamiltonian structure and conservation laws of the 

AMPLITUDE EQUATIONS 

Recalling the amplitude equations in the general form (|1.5p . i.e. 

M 

(A.l) idta„i = —i'd„i ■ V^-am + ^ ^ f^(p,q,r,m)0'pO-qO-rj 

S(Mp ,Mq - Air ) = S(Mm ) 

for TO = 1, ... , M, we want to highlight the Hamiltonian structure of this equation 
as well as to deduce the conserved quantities which follow from the specific structure 
of the resonances. For a general method on how to derive the reduced Hamiltonian 
structure and the first integrals for (|A.1[) we refer to [TTJ[T2]. Here we just collect 
the results that are obtained with the general theory developed there. 
We first note that (|A.1[) is a Hamiltonian system generated by 
„ M M ^ 

n''"^{a)^ Im(a„ i?™ • Vom) + -^^^^p^ apSqara™ dx, 

and the symplectic two-form i = Hence, (jA.ip takes the form 
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Note that the reduced Hamiltonian H"^'^{a) is not obtained as the lowest order 
expansion of the original Hamiltonian H^{u'^) of the full system, which reads 

Indeed, inserting the ansatz — w^, as given in (|4.ip . into TL°{u^) we find, as 
e — > 0, that 

„ M 

H'{u'j^)^I{a) + 0{e), with X(a) = / Ee^{k^)\a^\^ dx. 

•^R" m=l 

Even though not a Hamiltonian, X(a) clearly is a conserved quantity for the ampli- 
tude system ()A.ip . 

However, the energy levels Ei^{km) = ojm do not occur explicitly in the am- 
plitude system (jA.ip . Hence, any choice of ujm € R that is compatible with the 
resonance conditions (|2.4p leads to additional conserved quantities. More precisely, 
if lDi, . . . , ujm arc chosen such that for any p,q,r,m £ {1, . . . , M} the resonance 
conditions (j2.4p imply the identity 

COp — LOq -'r Ulr = ^m, 

then 

•^R' m=l 

defines a first integral. 

A particular choice is iii = • • • = 5^/ = 1, which leads to the trivial fact that the 
norm is preserved. The latter can of course also be deduced from the fact that 
the norm was preserved in the original problem for or from the fact that the 
mode system is invariant under the phase shifts (ai, . . . , um) ^ (e'"ai, . . . , e'"aAf), 
with a e M. Similarly, 1(a) can be understood as a first integral with respect to 
suitably chosen phase shifts that are compatible with the resonance structure, i.e. 

r„ : (ai,...,aA/)^ (e'=^^"ai,...,e'=^-"aA/),. 

Finally, it should be mentioned that the mode system (|A.1[) is translation invari- 
ant. This provides first integrals associated with the translation operators in the 
coordinate directions 

So far, we are not able to show that these conserved quantities are enough to 
provide a global existence result. Note in particular, that for the original problem 
the case k > leads to an energy that is definite and allows us to conclude global 
existence for the nonlinear Schrodinger equation p.ip . For the mode system (|A.1[) 
the sign of k is no longer helpful, since the term involving derivatives is indefinite. 
Note also that global existence for the mode system cannot be inferred from global 
existence of (jl.ip . since we can not expect the solutions to remain in the form 
of modulated pulses for all times. 
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